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FOREWORD 


This  report  was  prepared  by  Avco  Corporation/Systems 
Division  under  the  terms  of  the  Reentry  Antenna  Test  (RANT) 
Program,  Air  Force  Contract  No.  F04701  -68-C-G283.  It  is 
submitted  as  part  of  the  Preliminary  Design  Review  (PDR) 
package,  in  partial  fulfillment  of  the  contract  and  in  accordance 
with  Contract  Data  Requirements  List  (CDRL)  Item  No.  6.  The 
report  is  an  addendum  to  a  corresponding  report  submitted  at 
the  Systems  Studies  P.eview. 

The  program  is  administered  under  the  direction  of  the 
Space  and  Missile  Systems  Organization  (SAMSO)  with 
Captain  J.  Van  Zale  as  Project  Officer.  Mr.  C.  Okerstrom, 
Manager,  Aerospace  Corporation  Flight  Test  Vehicles  Project 
Office  and  Mr.  R.  A.  Mills,  Manager,  Aerospace  Corporation 
Electronic  Countermeasures  Project  Office  were  the  principle 
technical  monitors  for  the  work  reported. 
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1 . 0  Introduction 
-A 

'  In  Volume  2  of  the  RANT  preliminary  report  we  discussed  the  ground 
plane  model  for  predicting  the  effect  of  a  reentry  plasma  on  aperture  antennas. 
In  that  volume  we  also  began  the  discussion  of  the  cylinder  model,  (i.  e.,  aper¬ 
tures  on  a  plasma  covered  infinite  cylinder).  In  this  volume  we  conclude  the 
discussion  of  the  cylinder  model.  However,  for  completeness,  we  repeat  some 

i 

of  the  discussion  contained  in  Volume  2. 

In  modeling  the  plasma  the  effects  of  high  temperatures  have  been  inten¬ 


tionally  omitted.  Their  inclusion  would  substantially  increase  the  complexity 
of  the  calculations  without  significant  compensatory  gain  in  accuracy.  Analysis 
indicated  that  the  error  contribution  would  be  essentially  small  (only  a  few 
percentage  points)  at  the  electron  densities,  collision  frequencies,  and  tem¬ 


peratures  of  interest. 


In  both  models  the  plasma  is  approximated  by  a  dielectric  whose  local 
dielectric  constant  is  (all  analysis  and  results  assumes""1101  ,  except  for 
Tables  I,  II  and  HI  where  an  e  is  assumed  for  purposes  of  comparison) 


where 

€  =  local  dielectric  constant 

COp  =  electron  plasma  frequency 
Vc  =  electron-neutral  collision  frequency 
w  =  radian  signal  frequency 

The  collision  frequency  T^c  is  assumed  to  be  independent  of  field  strength. 
A  discussion  of  this  approximation  is  given  in  Section  3.0. 
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2 .  0  Cylinder  Linear  Plasma  Model 


2.  1  Analysis 


Tne  second  linear  plasma  model  considers  apertures  (either 
axial  or  circumferential )  cn  an  infinite  perfectly  conducting 
cylinder.  Again  we  assume  that  the  apertures  are  excited  in  the 
dominant  mode  and  that  the  inhomogeneous  plasma  is  approximated 
by  a  series  of  thin  layers,  each  of  which  is  individually  homogeneous, 
as  shown  in  Figure  l  a. . 

By  assuming  appropriate  expressions  for  the  fields  in  each 
layer  and  applying  the  boundary  conditions  on  the  tangential  fields 
across  the  boundaries  of  each  layer  and  across  the  aperture /neat 
snieid  interface,  we  can  obtain  the  following  expressions  for  the 
self  and  mutual  admittances 


CC 


A 

—  0 


/ 


•  C 

wWV 

/7k 


(2) 


where 


C  n 


(«) 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
fiwii  COPY  FUEMsHSB  TO  DDC 

-2- 


0.  O’  <4J 


IT  x> 


A 

C 


Af  ZA' 


**)  - 


-7 r  -os 


-  aperture  field  distribution  (dominant  mode) 


X,- 


-k  A  -  -*L  ?J. 


n 


n 


n 


n 


/ 

—  X. 

t  * 

/*• 

-±*J 

t  -  5 

/O*.  — Vn 

• 

. 

«/ 

^  -  •s . 

/>•  *''»v 

( 

V. 

(3) 

(4) 

(5) 

(6) 


<  I  r— »  - 

=  ^  ix  Xi-'v  -  r^)  \  (j±.*  j 
"4*44*0  K  ^44 
= fa  44 + (*)  h/44? 

^  4  hx  /4  *j  4  ^ 

=  Hn44^ 

-J-  X  4  [J/  (4  v~  Hn  X  v 

"4:  hJ4X  X4*) 

"  l_ -u  •  4^  ^  4  <4  4(^4  ^ J 

”  ^il  -4  (4  4 4  Mh  ,/ 

-  (/**/•)  Uk, 

=siv^ 4  HJ40A  (*) 

/A* 

-(*,*/*)  [x.  £4~)  (■<*;  HT/-44j 


£< 

a 


s# 

i* 

Of  I 

e-« 

|f 

M  |P 

!* 

TC 
A* 

co 

Si 


-5- 


THIS  PAGE  IS  BIST  QUALITT  PRACTICABJ^S 
tKvM  GO rX  K&KISiQS)  TO  BDG  _ . 


s 


SS 

ss 


B-i  Sj 

82 

3 

K?  S 

&  £ 

<0  if 


‘nm 


ex 


nm 


("filC'C+i'*'*. xxx  H JLJ 

-X^ 

—(v/j)  -^A>y-  ±L  ^ 

€ 


;  nm 


'•  /»*v  _  _  —j 

~  Hnfx*c)  X 

-  fir/z)  -^4^--^  jj^A  ^60  Hn^+4) 

X'V  — i 

'XXa  HXO  HnX  JJ 


¥ 


nm 


^  £0  /  '™rv 


* 


nm 


nm 


-U- 


£  £  y  /K-->n 


-  £=£— )  6 

*  Is  / 


/nw/ 


nm 


»*  *v 


~  )  2^- 


/r^V  0 


^m  -^<X  /3X- 


,  .  ./ 

m  +  1  -^L 


/O0V- 


/**•  ■*-! 

When  m  =  M,  which  corresponds  to  the  free  space  beyond  the 
plasma, 


«n  (M)  ~{fl 

n 


A 

0 

0 

ik*1 

0 

f!%M 

0 

■4^ 

0 

'fa 

0 

0 

-8- 


with 


nM 


/Atlf  OK 


dnM  “ 


nM 


f^nM 


-  (V/l)^  J-M  ^fA  \j^ b 0  Xm)  Hfl  ^/w) 

-^4  ^  '*'”)  ( X  m)J 

■£"£*.  X  CJL  ('X  J HJ Mo'*'") 

-  <«X£Vx„)xy^)] 


'IT  M 

2  £*/\ 


<3.  .  i 
f  i 


'i.4^  r 


nM 


X  HXxJlY^.X) 


-  £„-4  H„  fxM)  ££,'*'*)_ 

K  -t^  C4 2  >4  V  h  */xM)  x  fZ  ~*) 


~-^-u:-X)^)^(^ 


■nM 


h 


nM 


-  i  (^hJ(^  *»] 

-4  H»  (4**) £'(*«) 

-9- 


£t 

ns 


cy 


g  3 
«-<  ‘r 

S| 
8  3 

«  -3 
w  3 
** 
PH 


8 
£ 
W  : 


and  Hn(z)  =  Jn(z)  +  i  Y  n(z).  A  discussion  of  the  technique  used 

to  compute  the  Bessel  and  Neumann  functions  for  arbitrary  complex  arguments 
is  given  in  Section  2.3.  It  is  assumed  in  equations  (1)  and  (2)  that  the  apertures 
are  identical,  and  that  the  aperture  field  distribution  is  the  dominant  mode. 

For  the  case  of  two  axial  rectangular  slots,  as  shown  in  Figure  la,  with 
each  aperture  distribution  given  by 


(  IT*  \ 


equations  (1)  and  (2)  reduce  to 
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with  Kn  given  by  eqn  (6),  and 
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For  a  single  circumferential  gap  executed  by  a  uniform  z  -  directed 
electric  field,  as  shown  in  Figure  lb,  equation  (1)  becomes 
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where  Io  fcC)  is  given  by  equation  (3),  and  the  contour  of  integration  is  indented 
above  any  singularities  on  the  negative  Re(oC)  axis  and  below  any  singularities 
on  the  positive  Rc.(flC)  axis.  The  assumption  that  the  exciting  field  is  indepen¬ 
dent  of  0  ,  means  that  (9)  is  not  valid  for  cylinders  with  a  large  *A<L. 

The  single  axial  slot  radiation  problem  is  solved  by  applying  the  method 
of  steepest  descent.  Thus,  we  obtain  for  the  radial  component  of  the  Poynting 
vector 
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It  is  not  possible,  in  the  case  of  two  excited  axial  slots  on  the  cylinder, 
to  define  an  array  factor  as  was  done  in  the  ground  plane  model.  However, 
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the  radiation  pattern  (neglecting  mutual  coupling)  can  be  obtained  by  replacing 


the  term  e  {-tt,  °C)  in  the  definitions  of  Ln  and  Nn,  above,  by  the  term 

A  *  ,  T  _  l  C'Xt»  +  m  f .  4  »<  2|)1 


where 


^  =  phase  difference  between  the  excitations  of  the  two  slots 

For  the  circumferential  gap  the  radiation  pattern  is  given  by 


#  1 
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Equations  (10)  and  (11)  are  not  valid  for  0  =  o.  This  is  because  in  deriving 

(10)  and  (11)  we  have  used  the  asymptotic  expansion  for  Hn(-^^tsin  0)  which 

is  not  valid  for  0  =  o  .  (In  addition  as  0^o  the  saddle  point  approaches  a 

singularity,  so  that  the  usual  saddle  point  expression  would  not  be  valid  in 

\ 

any  case.  ) 

2.2  Discussion  of  Cylinder  Linear  Plasma  Compute,r  Program 

Using  all  the  above  results,  Avco  Computer  Program  2817  was  developed. 

It  has  four  options  as  follows: 

(a)  Option  1  computes  ,  Yj^,  and  isolation  for  a  pair  of  identical 
axial  slots  on  an  infinite  cylinder  covered  by  an  arbitrary  inhom- 
geneous  plasma.  (The  plasma  may  be  lossy  or  lossless. ) 

(b)  Option  2  computes  the  radiation  pattern  for  both  a  single  and 
axial  slot  and  for  a  pair  of  axial  slots. 
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(c)  Option  3  uses  eqn  (9)  to  compute  the  admittance  of  a  circumferential 
gap  on  an  infinite  cylinder  covered  by  an  arbitrary  inhomogeneous 
plasma. 

(d)  Option  4  calculates  the  radiation  pattern  for  the  above  situation. 

A  complete  Fortran  Listing  for  program  2817  is  contained  in  a  separate 
volume.  A  typical  input  sheet  and  explanation  of  inputs  is  contained  in 
Appendix  III. 

2.2.1  Verification  and  Results  from  Options  1  and  2 

To  verify  the  results  of  program  2817,  we  have  compared  its  results  with 
some  results  obtained  by  NASA,  and  also,  where  applicable,  with  ground  plane 
results. 

Table  I  compares  the  values  provided  by  a  computer  program  with  similar 
values  provided  by  a  computer  programs  developed  by  C.  Knop  (1965)  for  the 
NASA  Langley  Research  Center.  The  antenna  considered  in  this  case  was  a 
.ixial  rect.i.igular  slot  (<x  0O  =  .  1687  and  b  =  2.  38)  located  on  an  infinite  cylinder 
of  radius  21  =  1  covered  by  a  homogeneous  plasma  of  thickness  ,  and 
dielectric  constant  £  One  may  note  from  Table  I  that  the  comparison  is 
reasonably  good.  The  study  presently  being  performed  by  Golden  and  Stewart 
(1969)  of  Aerospace  will  afford  additional  test  cases  in  the  future. 

To  examine  the  case  when  6 is  complex  we  take  A  =  3.  6,  b  =  8.  5,  k  =  .418, 
^  I  =  1  +  i  20,  and  A  ^  =  .  0534.  In  this  case  we  expect  that  the  thin  sheath 
approximation  is  valid  (Fante,  1967), and  that  Yjj  =  YjJ0)  +  Ys,  where  Yj^°^ 
is  the  free  space  admittance  and  Yg  =  i  kA^  (€  -  1)(1207T )”  *.  (The  error  in 
using  the  above  formula  for  the  present  example  is  of  order  12%).  Using 
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TABLE  I 


COMPARISON  OF  ADMITTANCES  FOR  AN  AXIAL 
SLOT  ON  A  PLASMA  COVERED  CYLINDER 


*1 

j  CONDUCTANCE 

!  k  Avco^  NASA^ 

j  SUSC.EPTANCE 
Avco1  NASA2 

.01 

*  / 

'  1.5  !  1.  50  X  10"4 

i 

1. 48x  10-4 

1 

-1.  54xl0-4 

-1.  52xl0-4 

.01 

1.  5  1.  52  x  10"4 

t 

1.52x10-4  } 

-1.  92xl0"4 

-1.  88x10-4 

From  Avcg  computer  program  2817 

NASA  Langley  Computer  Program  fay  C.  Knop 


the  thin  sheath  approximation  we  find  that 

Y  =  1.33  x  1 0  "  3  -  i .  2 10  x  1 0“  3  mho 

while  program  2817  calculates 

Y  =  1.27  x  10“ 3  -  i  .24  x  10*3  mho 

Since  the  value  calculated  by  2817  agrees  with  the  thin  sheath  result  (to  within 
the  error  inherent  in  the  thin  sheath  approximation)  we  conclude  that  program 
2817  is  also  operating  satisfactorily  for  the  case  of  complex  £  . 

It  is  also  desirable  to  study  the  isolation  results  predicted  by  program 
2817.  Figure  2  shows  the  isolation  of  two  axial  slots.  The  upper  curve 
shows  the  change  in  H-plane  isolation  as  a  function  of  separation,  while  the 
lower  curve  describes  the  change  in  E-plane  isolation  as  the  azimuthal 
separation  is  varied.  As  expected  the  isolation  increases  with  increasing 
slot  separation. 

Figure  3  shows  the  cnange  in  E-plane  isolation  for  a  cylinder  covered 

by  a  two  layer  plasma.  For  this  case  the  computed  self-admittance  is 

Yu  =  .  758  x  ID'2  +  i  .  118x  10-2.  It  should  be  pointed  out  that  in  this 

example  the  sum  over  n  (the  azimuthal  modes)  was  terminated  at  N  =  101 

due  to  computer  underflows.  Generally,  for  small  isolations  (i.  e. ,  isolations 

less  than  about  80  db)  this  does  not  present  any  problem.  However,  for  cases 

when  the  isolation  is  large,  and  hence  Y ^  very  small,  it  is  possible  that  the 

© 

....  )  may  be  of  the  same  order  as  the 

to 

computed  answer.  Then,  the  computed  isolations  should  be  regarded  in  an 


neglected  terms  (i.e.,  21 
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order  of  magnitude  sense,  and  may  possibly  be  a  few  ab  different  from  the 

answer  one  would  obtain  if  the  T  were  not  truncated  t  '01.  {The  neglected 

TV 

terms  are  especially  important  for  cases  when  0,  is  largo  since  there  then 

arises  cancellation  for  small  n  from  term  to  term  in  the  su  ->at:on 

el«01  F  (  m)  in  eqn  (8),  so  that  the  large  n  terms  assume  increasing 
n 

importance. )  Fortunately,  in  practical  situations  we  do  not  need  extremely 
accurate  results  when  the  isolation  is  large. 

It  is  also  interesting  to  conjecture  how  the  H-plane  admittance  and 
isolation  (for  this  case  01  =  o)  between  two  axial  slots  separated  by  a  distance 
Zj  on  a  cylinder  compare  with  the  admittance  and  isolation  between  the  same 
slots  (separated  by  the  same  distance  zj)  on  a  ground  plane.  A  comparison  is 
shown  in  Table  II.  We  note  that  the  self -admittance  Yjj  for  the  slot  on  the 
cylinder  does  not  differ  very  vastly  from  that  for  the  same  slot  in  a  ground 
plane,  even  though  the  cylinder  circumference  is  only  1.  5  wavelengths. 

Even  the  difference  between  the  values  of  isolation  is  not  very  large. 

Finally,  let  us  consider  the  radiation  pattern.  Figures  4a  and  4b  show- 
some  typical  9  cuts  (for  0  =  oc)  for  both  a  single  slot  and  a  pair  of  H-plane 
slots  which  are  excited  out  of  phase,  but  w'ith  equal  amplitudes.  Figure  4a 
shows  the  case  when  the  cylinder  is  covered  by  a  layer  of  lossy  plasma 
and  a  dielectric.  Figure  4b  shows  the  case  when  the  cylinder  is  covered  by 
a  very  slightly  lossy  dielectric  layer  separated  by  a  thin  air  gap  from  the 
cylinder  surface.  Note  that  in  this  case  the  pattern  is  peaked  near  9  =  0, 
indicating  the  effect  of  surface  wave  modes. 


comparison  of  Results  for  Slots  on  a  Ground  Plane 
with  those  on  a  Cylinder  for  a  =  3.6,  b  =  8.  5, 


2.2.2  Verification  and  Results  from  Options  3  and  4 

To  verify  the  validity  of  c  ptions  3  and  4  we  have  first  considered  the 
limiting  case  of  a  bare  cylinder  (for  geometry  see  Figure  lb).  For  this 
problem  some  values  have  been  computed  by  Swift  et  al  (1969)  of  NASA. 

Table  III  is  a  comparison  of  their  results  with  those  from  option  3  of  program 
2817.  We  note  that  the  agreement  is  excellent.  Also  for  ka  «  1  we  can 
compare  2817  with  the  theoretical  admittance  formula  of  Fante  (1966)  and 
Miller  (1967)  given  by 


(12) 

where  YQ  is  the  admittance  of  free  space  and  T  =  1.781.  The  results  of 
2817  have  been  compared  with  (12),  and  the  agreement  was  excellent. 

To  verify  the  validity  of  option  3  when  the  cylinder  is  covered  by  a 
plasma,  we  have  again  chosen  a  plasma  which  meets  the  conditions  necessary 
for  the  thin  sheath  approximation  to  be  valid.  We  take  a  =  3,  k=  1,  b=.  105, 

A  i  =.02  and  6  j  =  1  +  i  2 5 .  The  input  admittance  as  predicted  using  the 
thin  sheath  approximation  is 

Yn  =  1.49  x  10  "3  -  i  3.44  x  10  "4  vrAo 

while  the  value  obtained  from  program  2817  is 

Yjj  =  1.  35  x  10"3  -  i  4.0  x  10'4  *mho 

Again  the  agreement  is  favorable,  and  the  differences  are  within  those 
expected  due  to  errors  in  using  the  sheath  approximation. 
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fun 


Comparison  of  admittances  for  an  axial  slot 
on  a  bare  cylinder 


Yj  ^  K* ) 


kJ 

Limit  of. 28 17 

1 

.  035 

.685  x  10"4  +  j  1.49  x 

10'4 

2 

i 

.070  * 

i 

^  1.  16  x  10-4  +  j  2.  52  x 

10-4 

Yjj  (vrU) 

NASA 

.685  x  10"4  +  j  1.49  x  10-4 
1.  16  x  10"4  +  j  2.  54  x  10-4 


I 

$ 

3 

i 

| 


3  .  105  .  1.62  x  10"4  +  j  3.43  x  10  “4  1 . 62  x  10‘4  +  j  3. 44  x  10"4 

,  i  w 


To  verify  the  radiation  patterns  as  predicted  by  option  4,  we  again 


first  consider  the  free  space  limit.  In  this  case  (11)  reduces  to: 

\  \  2 


R  = 


1 


/  /  kb  A  ^ 

/  sin  (  <T  cos  6  j 


sin‘  0  fHo  (ka  sin  0)1  ^  j  kb 

I-  J  '■  —=r~  cos  9 


/ 


/ 


(13) 


A  comparison  with  the  results  of  (13)  with  the  results  computed  by  option  4 
of  2817  is  given  in  Table  IV. 

To  verify  the  validity  of  the  computed  radiation  pattern  when  a  plasma 
is  present,  we  consider  a  plasma  for  which  the  thin  sheath  approximation  is 
applicable.  In  that  limit  it  can  be  theoretically  shown  (Fante,  1967)  that  if 
the  aperture  excitation  is  unchanged  the  pattern  is  also  unchanged  (both  in 
amplitude  and  shape,  except  near  0  =  o)  from  the  vacuum  situation  when  a 
plasma  layer  is  present.  Table  V  shows  the  ratio,  R^,  of  the  radiation 
pattern  with  a  plasma  (satisfying  the  thin  sheath  criterion)  to  that  in  vacuum. 
We  note  that  the  ratio  is  generally  near  unity  as  is  expected  from  the  above 
theoretical  considerations. 
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TABLE  IV 


Comparison  of  Theoretical  and  Computed  Vacuum  Radiation  Patterns  for 
k  =  .25,  a  =  4,  b  =  .  140.  Values  are  Normalized  to  Unity  at  0  =90°. 


Q  (degrees) 


R  from  2817 


R  from  equation  (13 


TABLE  V 


Ratio  Rp  of  the  Radiation  Pattern  with  a  Plasma  Layer  (with  €,  =  1  +  i  25 
A  i  =  0.  1)  to  the  Vacuum  Pattern  for  a  =  4,  b  =  .  140,  k  =  .25 


j^p_ 

6  (deg 

.  999 

90 

.  998 

80 

.  996 

70 

.991 

60 

.  982 

50 

.  967 

40 

.  936 

30 

.  860 

20 
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2. 3  Techniques  for  Computing  the  Bessel  and  Neumann  Functions 


Our  method  for  computing  Bessel  functions  is  based  on  a  technique 

described  in  the  S.I.  A.M.  Review,  Vol.  9  No.  1  January  1967,  by 

Walter  Gauschi.  The  article,  entitled  "Computational  aspects  of  three- 

term  recurrence  relations"  contains  a  complete  description  of  the  method 

and  gives  a  full  error  analysis.  The  program  will  produce  highly  accurate 

Bessel  functions  provided  the  imaginary  part  of  the  argument  is  less  than 

174.673  in  magnitude,  and  therefore  no  asymptotic  expansions  are  necessary. 

Limiting  forms  are  used  for  arguments  whose  magnitude  is  less  than 
-8 

1.  x  10  ,  although  this  is  not  essential  but  merely  faster. 

The  Neumann  functions  are  computed  using  recursion  upward 
(a  numerically  stable  procedure)  with  starting  values  for  Y0(z)  and  Yj(z) 
obtained  by  Neumann's  expansion  for  Yn  in  terms  of  Jj ,  i  =  1,  2,  . . . 

Hankel  functions  H^^(z)  are  computed  as  J*(z)  +  i  Y^(z)  unless 
a  loss  of  significance  of  six  or  more  places  will  result.  This  occurs  only 
for  large,  pure  imaginary  arguments.  Here  asymptotic  expansions  for 
H^^(z)  and  (z)  are  used.  A  reference  for  the  Y  and  H  expansions  is 

the  Handbook  of  Mathematical  Functions  published  by  the  National  Bureau 
of  Standards,  U.S.  Dept,  of  Commerce. 


. . . . . . 


3.0  Discussion  of  the  Collision  Freauencv  Approximation 


In  programs  2656  and  2817  we  have  neglected  the  effect  of  electric 


field  strength  on  the  electron-neutral  collision  frequency.  It  is  readily 


shown  that  for  the  range  of  altitudes  over  which  the  signal  attenuation  is 


significant,  this  approximation  is  appropriate.  Figure  5  shows  the 


effective  collision  frequency,  calculated  from  the  Boltzmann  equation  model, 


as  a  function  of  E  /p.  For  the  2  watt  RANT  S-Band  slot,  in  the  altitude 


regime  between  130  kft  and  40  kft,  the  slot  conductance  varies  between 


10'^  and  2  x  10“^  mhos/meter.  As  a  result  the  electric  field  outside  the 


antenna  is,  at  most,  of  order 


1/2  [ 


(2)  (2) 


|  1  -  *  — * 

-  (2.  54)(4.45)  x  10-3 


18.  7  v/ cm 


At  130  kft  the  pressure,  assuming  T  =  4000CK.  is  of  order  30  torr.  For 


lo-wer  altitudes  the  pressure  will  be  higher.  Thus,  in  the  regime  from 


to  to  130  kft,  where  attenuation  is  most  significant  (and  an  accurate  deter¬ 


mination  of  important)  the  maximum  value  of  Eg  /p  is  of  order 


Referring  to  Figure  A.  1  we  see  that,  since  /p  <  1/2  for  the  2  watt 


communications  link,  there  is  almost  no  error  in  approximating 


(E)  by  v'c  (E  =  0). 
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Program  2656  -  Linear  Ground  Plane  Plasma  Model 


Program  2656  calculates  the  admittance,  isolation  and  radiation 
pattern  for  rectangular  slots  in  a  ground  plane  covered  by  an 
outwardly  inhomogenous  plasma. 


Program  2656 
EXPLANATION  OF  INPUT 


NUMD 

NLONG 

NOP 


ER 


A 

B 

XAMBDA 

RAD 

D 


Number  of  dielectric  slabs.  Maximum  value  is  100. 
Must  be  >  1  if  NOP  is  1  or  3  or  4. 

Optional  input.  NUMD  may  be  trucated  to  conserve 
computer  time  in  admittance  calculation.  NLONG  is 
then  the  total  number  of  dielectric  slabs. 


Option: 


if  =  1,  then  integral  for  admittance  with  plasma 

(Y*  ). 

o 

if  =  2,  then  integral  for  admittance  with  free 

space  ( Y*  ). 

ov 

if  =  3,  both  option  1  and  2  with  complete  printout 
of  radiation  pattern. 


if  =  4,  radiation  pattern  with  admittance  calculation: 
can  include  admittance  by  inputting  XKP  and/or  XKV, 

if  =  5,  expressions  for  near  fields  of  a  slot 
antenna  covered  by  a  heat  shield.  Need  x,  y 
and  z.  If  E^,  E  ,  or  E^  =  0  then  X0P,  Y0P 
or  Z0P  must  be^'O  respectively. 


Used  in  formula  as  follows: 

T 


Yc  = 


t 


ER  -  ( 


a  k 


>■ 


r 


/  377 


ER  is  relative  dielectric  constant  of  material  inside  wave 
guide.  Preset  to  1. 

Large  wave  guide  dimension  in  Cm. 

Small  wav  .•  guide  dimension  in  Cm. 


Signal  wave  length  in  centimeters. 

Factor  defined  D  =  D»RAD 

Thickness  of  each  dielectric  layer  in  Cm.  Maximum  of 
100  values. 
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EPSR  Relative  heat  shield  and  plasma  dielectric  constant 

(real  part).  Maximum  of  100  values.  Not  used  in 
N0P2. 

EPSI  Imaginary  part  of  dielectric  constant.  Vector  as  above. 

Not  used  in  N0P2. 

NOR  Number  of  R's.  IS  NOR  <  100;  used  in  N0P  3  or  4. 

NTHETA*  Number  of  thetas,  1  <  NTHETA  <  360;  used  in  NOP  3  or  4. 

NPHI  Number  of  phis.  1  <  NPHI  <  360;  used  in  NOP  3  or  4 

R  Vector.  Distance  from  antenna  slot  in  meters.  Maximum 

of  100  values.  Used  for  N0P  -  3  or  4  (preset  to  1.  ) 

THETA  Vector.  Polar  angle  in  degrees.  (zero  not  permitted) 

Maximum  of  360  values;  used  for  NOP  3  or  4. 

0  <  THETA  <180;  theta  *  90°. 

PHI  Vector.  Azimuth  angle  in  degrees.  (zero  not  permitted) 

Maximum  of  360  values;  used  for  NOP  3  or  4, 

0  <  PHI  <  360;  No  multiples  of  90°. 

XKP  Multiplying  factor.  |l  +  Rp[^  times  XKP  replaces 

1 1  +  Rp|2  in  NOP  3  or  4. 

However,  in  N0P4  Rp  =  0.  (XKP  preset  to  1.) 

XKV  Multiplying  factor  |l  +  Ry|^  times  XKV  replaces  |l  +  RyJ^ 

in  N0P  3  or  4.  However,  in  NOP  4  Ry  =  0. 

(XKV  preset  to  1. ) 

IREND  Number  of  successively  larger  R's  to  be  used  in  integrals 

(Preset  to  20). 

RINFCN  Convergence  criterion  for  integrals  as  R-**90  ,  preset  to  .01 

NPLT2P  Plot  option  for  (N^)  2p  . 

If  =  1  plots;  if  =  0  no  plots,  (preset  to  0) 

NPLTNV  Plot  option  for  (N  ),r. 

r  V 

If  =  1  plots;  if  =  0  no  plots.  (Preset  to  0) 

Plot  option  for  (N  ),  /  (N  ),.. 

r  &p  r  V 

If  =  1  plots;  if  =  0  no  plots.  (Preset  to  0) 
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NPLTN2 


l6gr 

NV2* 

TITNRP 

TITNV 

TITN2 

EPSIH 

EPSRH 

PG 

R A  TEMP 
RADEN 
EL 
GL 

PL 

NL 

TEMP 

DENSRA 

EMC 

ELEDEN 

C0LFRE 

El 


Option:  if  =  1  then  polar  plot  of  R  is  linear, 
if  =  2  then  log^R,  (preset  to  1). 

Option:  if  =  1  then  points  of  plots  are  not  connected, 
if  -  2  then  connected  (preset  to  2). 

56  character  title  for  (Nr)£p  polar  plot.  The  title  must 
appear  on  the  second  card.  Requires  2  cards  if  used. 

56  character  title  for  (Nr)y  polar  plot.  Requires  2  cards  if  used. 

56  character  title  for  (Nf)2p/  (RTr)y  P°*ar  plQt*  Requires  2  cards 
if  used. 

Optional  input.  Imaginary  part  of  heat  shield  dielectric 
constant.  Can  be  used  only  if  IN0P  f  3. 

Optional  input.  Real  part  of  heat  shield  dielectric  constant. 

Boundary  layer  pressure,  atmospheres. 

Factor  multiplying  temperature  points.  Preset  to  1. 

Factor  multiplying  density  ratios.  Preset  to  1, 

Ionization  potential ,  J[  th  electrophillic  negative  ion. 

Ratio  of  electronic  weighting  factors  of  each  electrophillic 
species. 

Partial  pressure  of  Jt  ^  electrophillic  species. 

Number  of  electrophillic  species.  Maximum  of  10. 

Temperature  of  each  plasma  layer,  deg  K,  maximum  of  99  values. 
Density  ratio  of  each  plasma  layer. 

Ratio  of  heat  shield  density  to  total  gas  density  in  each  plasma  layer. 
Electric  density  of  each  plasma  layer,  particle /cm  . 

Collision  frequency  of  each  plasma  layer,  rad / sec. 

Ionization  potential  of  i  ^  contaminant,  ev. 
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EMI 


Parts  Per  Million  of  i  ^  contaminant. 


GI  ^ .  ratio  of  electronic  weighting  factors  of  each 

gi 

contaminant. 

WI  Atomic  weight  of  each  contaminant 

NI  Number  of  contaminants.  Maximum  of  10. 

IN0P  3  input  EPSR  and  EPSI 

2  input  ELEDEN  and  C0LFRE 
1  input  ELEDEN,  TEMP  and  DENSRA 
0  input  TEMP  and  DENSRA  (program  computes  clear  air 
values  of  ELEDEN  and  CDLFRE) 

-1  input  TEMP,  DENSRA,  EMC,  EMI,  El,  GI,  WI,  NI 
(program  computes  ELEDEN  only  due  to  contaminants 
and  COLFRE  for  air) 

-2  input  same  as  -1  (program  includes  effect  of  air  on  ELEDEN) 

IS0L0P  Option:  If  set  to  1  then  the  factor  cos  (A  sin  8)  cos  (B  cos  8) 

is  also  used  in  the  integrand  for  N0P  =  1  and/or  N0P  =  2P 

AO  The  H  plane  slot  separation  in  cm.  Needed  for  IS0L0P  =  1 

and  GAMMA  calculation  in  N0P  3  and  4.  Preset  to  0. 

BO  The  E  plane  slot  separation  in  cm.  Needed  for  ISDL0P  =  1 

and  GAMMA  calculation  in  N0P  3  and  4.  Preset  to  0* 

EPSIL  When  Z0P  =  1»  and  EPSR  =  1*  and  EPSI  =  0*  ,  the  relative 

interval  (  (  )  in/»  used  at  singularity.  Also  used  for 
N0P  2.  Preset  to  .  1 

N<t)XYZ  The  number  of  x,  y  and  z's  to  be  inputted. 

X  t  Vectors.  In  E  (x,  y,  z),  E  (x,  y,  z)  and  E2(x,  y,  z). 

Y  ,  Used  only  in  NDP  =5  Y 

Z  j  Maximum  of  100  values 

X(t)P  Option.  If  =  0»  then  convergence  not  used  in  solution  of  Ex, 

If  =  1*  then  convergence  used.  Preset  to  1*  If  Ejj  =  0 
X(J)P  must  be  0. 


A-8 


Y0P 

Z0P 


Option.  Same  as  X0P  but  explanations  apply  to  E 
and  E^  respectively.  ^ 

N0DBLE  The  maximum  number  of  doubling  of  points  in  the 
Simpsons'  integral  solution.  Preset  to  8. 

If  greater  accuracy  required,  use  10  or  12  but  problem 
will  increase  in  time  radically. 

THTCNV  Theta  convergence  criterion.  Preset  to  .  005 

IS  KIP  Option:  If  =  1  and  if  IS0L0P  =  1  this  causes  the  suppression 

of  the  Y  calculation.  Y  value  used  in  the  IS0L  calculation 
can  either  be  inputted  or° the  last  value  from  a  previous  case 
will  be  used.  (Preset  to  0) 

YREAL1  The  real  part  of  Y  for  N0F  =  1.  Used  when  ISKIP  -  1 

o 

YIMAG1  The  imaginary  part  of  Y  for  N0P  =  1.  Used  when  ISKIP  =  1 

Y REAL 2  The  real  part  of  Y  for  N0P  =  2.  Used  when  ISKIP  =  1 

o 

YIMAG2  The  imaginary  part  of  Y  for  N0P  =  2.  Used  when  ISKIP  =  1 

o 

Note.  Calculated  values  of  ELEDEIv  and  C0LFRE  assume  equilibrium, 
EPSHR  and  EPSIH  may  be  used  only  when  IN0P  ^  3. 
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OUTPUT 


The  input  program  outputs  the  card  images  of  the  input. 


OUTPUT  NAME 


DEFINITION 


YO'  Real  and  Imaginary  part  of  Self  Admittance 

Y12  Real  and  Imaginary  part  of  Mutual  Admittance 

IS0L  Isolation  in  DB 

YG  Wave  Guide  Characteristic  Admittance 

YOV  Vacuum  Self  Admittance 

(NR)P  Radiation  Pattern  with  Plasma 

(NR)V  Vacuum  Radiation  Pattern 
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OPERATION  INSTRUCTIONS  FOR 
PROGRAM  2656 


Program  26-A  computes  the  self*  and  mutual  admittances  of  unmatched 
rectangular  slots  on  a  plasma  covered  ground  plane.  Attached  here  is 
an  input  sheet  and  a  complete  explanation  of  all  inputs.  We  have  also  inclu¬ 
ded  three  sample  outputs.  Exhibit  A  is  a  case  where  ti  e  admittance  and 
isolation  are  computec  for  a  two  layer  lossy  dielectric  (or  plasma).  Ex¬ 
hibit  B  computes  the  admittance  and  isolation  for  a  ground  plane  covered 
b\  a  slightly  lossy  dielectric.  Vv  note  that  in  this  case  the  program  prints 
out  a  statement  that  the  integral  c.d  not  converge  in  a  given  region  duo  to 
the  tight  convergence  criteria  we  have  imposed.  If  the  real  parts  are  within 
about  of  each  other  and  the  imaginary  parts  are  also  within  about  5°)  of 
each  other,  the  error  in  Yu  and  Y]2  is  usually  less  than  1%.  If  the  real 
or  imaginary  parts  differ  by  more  than  5%,  the  user  should  increase 
NG5DBLE.  Exhibit  C  illustrates  the  case  w'hen  the  admittance,  radiation 
pattern  and  attenuation  are  calculated  for  a  multilayer  plasma. 

Also  included  in  this  package  is  the  complete  deck.  This  deck  is  appropr ,  -tu¬ 
tor  use  on  the  IBM  360-75. 


*Note  that  e-lw‘  dependence  is  assumed.  Thus  2656  values  for 
admittance  are  complex  conjugate  cf  those  for  e+iwt. 
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PURPOSE  OF  PROGRAM  28l? 


Cylinder  Linear  Plasma  Model 


Program  2817  confutes  the  admittance  isolation  and  radiation  pattern  for 
slots  on  an  infinite  plasma  covered  cylinder.  The  program  has  four  options: 
Option  1  calculates  the  admittance  and  isolation  for  a  pair  of  axial  slots, 
option  2  calculates  the  radiation  pattern  for  the  above  case,  option  3  cal¬ 
culates  the  admittance  cf  a  circumferential  gap,  and  option  i  calculates 
the  radiation  pattern  for  the  circumferential  gap. 
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DIGITAL  COMPUTER  IHPUT 

PROBLEM  NO. 

MEMO  NO. 

SECTION  NO. 

CONTINUATION  SHEET 

REQUEST  FORM 

281? 

EXPLANATION  OF  INPUTS 

A  Radius  of  the  bare  cylinder. 

B  Slot  width  in  the  axial  direction. 

BETA  Amplitude  of  the  excitation  of  slot  2  relative  to 

that  of  slot  1, 

COMPARE  When  running  option  4,  shall  the  pattern  be  compared  to 
one  for  £  =  {1,  0)  .  Answer  .  YES  or  .  N0 

DELTA1  Thickness  of  the  heat  shield. 

DELTA  Incremental  thickness  of  the  plasma  layers. 

EPSIL0N  Relative  dielectric  constant.  (List  the  real  and  imaginary 
part  for  each  layer  covering  the  bare  cylinder. ) 

K  Free  space  wavenumber. 

LAMBDA  Phase  of  slot  2  relative  to  slot  1. 

PHIO  Angular  width  of  the  aperture  in  degrees. 

PHU  Azimuthal  separation  of  the  centers  cf  slots  1  and  2, 

in  degrees. 

Z1  Axial  separation  of  the  centers  of  slots  1  and  2, 

in  degrees. 

^/LAYERS  Total  number  of  layers  (heat  shield  plus  plasma). 

PHI  Value  or  values  of  ^  for  radiation  pattern. 

THETA  Value  or  values  of  &  for  radiation  pattern.  Values 

may  be  listed  or  given  in 

STAR  TIN G(INCREMENT)FINAL 

format.  EXAMPLE:  1.  (1.  )  10.  (2.  )80,  (1.  )89. 

/>PTI$N  Use  1  if  we  desire  a  calculation  of  YT1,  Y12  and  isolation 
for  axial  slot. 

Use  2  if  we  desire  radiation  pattern  only  for  a  single  slot 
and  a  pair  or  axial  slots. 

Use  i  if  we  desire  a  calculation  of  Yll,  Y12  and  isolation 
for  a  circumferential  slot. 


Use  4  if  we  desire  the  radiation  pattern  of  a  circumferential 
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EXPLANATION  OF  OUTPUT  SYMBOLS 


The  self  admittance 


ine  mutual  admittance 


The  waveguide  characteristic  admittance 


The  isolation  in  db 


The  far  radial  Povnting  vector  of  a  single  slot 


The  far  radial  Poynting  vector  for  a  pair  of  slot, 
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PURPOSE  OF  PROGRAM  OK- 5 


Program  OK- 5  was  written  to  modify  the  attenuation  and  isolation  calculated 
by  Avco  computer  program  2656  to  account  for: 

1.  The  fact  that  the  S-  and  C-band  experiment  antennas  are  not 
infinite  waveguides  but  cavity  backed  slots. 

2.  The  fact  that  there  may  be  a  lossless  tuning  network  in  the 
transmission  line  feeding  the  probe. 

Multiple  cases  can  be  run  by  stacking  one  case  behind  the  other.  Ho 
"//  XEQ"  card  is  necessary  between  cases.  After  the  last  data  card  of 
the  last  case  being  run  a  card  with  a  "1"  punched  in  column  two  should  be 
used  in  order  to  properly  terminate  program  execution. 


C-3 
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Flow  Chart  for  Program  OK- 5 


Input  Fbrmat  Description  For  OK- 5 


Input  Description 


Code  Word  When  set  to  1 
program  stops 

Title  of  Case 

Attenuation  outputed  from 
2656  when  reference 
antenna  is  radiating 
into  the  reference 
environment 

|l+RP|  **2  Transmission 
Coefficient  for  A2656 
Attenuation 

Attenuation  Outputed  From 
2656  When  Subject 
Antenna  is  Radiating 
Into  Present  Environ¬ 
ment 

jl+RP|  **2  Transmission 
Coefficient  for  B2 656 
Attenuation 

Complex  Aperture  Self-Ad¬ 
mittance  to  Which  Re¬ 
ference  Antenna  is 
Matched 

Complex  Aperture  Mutual- 
Admittance  to  Which 
Reference  Antenna  is 
Matched 

Complex  Aperture  Self -Ad¬ 
mittance  into  Which 
Reference  Antenna  is 
Radiating 

Complex  Aperture  Mutual- 
Admittance  into  Wnich 
Reference  Antenna  is 
Radiating 


IT 

A2656  2 


Y120R 


Columns 

Format 

Typical 

Value 

1-2 

12 

0 

3  - 

20A2 

TEFLON,  etc. 

1-10 

F10.0 

-0.7369 

11  -  20 

F10.0 

2.28611 

1-10 

F10.0 

-0.701*3 

11  -  20 

F10.0 

2.27707 

1-20 

2F10.0 

0.001071 

O.OOO5976 

21  -  4o 

2F10.0 

0.0 

0.0 

1-20 

2F10.0 

0.001071 

0.0005976 

21  -  1*0 

2F10.0 

0.00004792 

-0.0001681 

Input  Description 

Symbol 

Card 

Columns 

Format 

Typical 

Value 

Complex  Aperture  Self-Ad¬ 
mittance  to  Which 

Subject  Antenne  is 

Matched 

Y110 

6 

1  -  20 

2F10.0 

0.001071 

O.OOO5976 

Complex  Aperture  Mutual- 
Admittance  to  Which 
Subject  Antenna  is 

Matched 

Y120 

6 

21  -  1*0 

2F10.0 

0.0 

0.0 

Complex  Aperture  Self-Ad¬ 
mittance  into  which 
subject  antenna  is 
Radiating 

Yll 

7 

1-20 

2F10.0 

0.001081 

0.0005911 

Complex  Aperture  Mutual- 
Admittance  into  which 
Subject  Antenna  is 
Radiating 

Y12 

7 

21  -  1*0 

2F10.0 

0.0000171*5 

-0.0001671* 

Aperture  Length  (cm) 

A 

8 

1  -  10 

F10.0 

1*.1*1* 

Aperture  Width  (cm) 

B 

8 

11  -  20 

F10.0 

2. 51* 

Length  of  Probe  in  Cavity 
(cm) 

D 

8 

21  -  30 

F10.0 

2.0 

Probe  Radius  (cm) 

R 

8 

31  -  1*0 

F10.0 

0.159 

Distance  From  Probe  to  Back 
of  Cavity  (cm) 

ALL 

8 

1*1  -  50 

F10.0 

2.0 

Distance  From  Probe  to  Aper¬ 
ture  (cm) 

ABL 

3 

51  -  60 

F10.0 

2.1*1* 

Relative  Dielectric  Coef¬ 
ficient  for  Cavity 

EPS 

9 

1-10 

F10.0 

3.75  . 

Characteristic  Admittance  of 
Line  Feeding  the  Probe 

YP 

9 

11  -  20 

F10.0 

0.02 

Free  Space  Wavelength  of 
Signal  (cm) 

AIMDA 

9 

21  -  30 

F10.0 

12.766 

OK- 5  Output  Description 


Description 

Input 

Magnitude  of  the  reflection  coefficient  squared  for 
the  reference  antenna. 

Magnitude  of  the  transmission  coefficient  squared 
for  the  reference  antenna. 

Reflection  loss  in  dB's  for  the  reference  antenna. 


Isolation  in  dB's  for  the  reference  antenna. 


Input  y 

Magnitude  of  the  reflection  coefficient  sqared 
for  the  subject  antenna. 

Magnitude  of  the  transmission  coefficient  squared 
for  the  subject  antenna. 

Reflection  loss  in  dB's  for  the  subject  antenna. 

Isolation  in  dB’s  for  the  subject  antenna. 


Attenuation  loss  in  a3’s.  Given  by 


Attenuation  =  -IOLOG-^q 


/Power  radiated  by  subject  antenna  v 
'Power  radiated  by  reference  antc;.na/ 
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OK- 5  OUTPUT 

®  TEFLON  S-RAND  150KFT  NEAR-H 


'SLOT  PARAMETERS 


SLOT  LENGTH  =  4.440 

PR0°F  DEPTH  =  2'*D"00" 

D~0=5E  TO  B A C<  =  2.000 

-poC°E  ’ADVITT .  =  UTOTiCT 

WAVELENGTH  =  12. ^66 


slot  WIDTH  : 
"PROBE  PAD  I  US — 5 
PROBE  TO  FRONT 


2.540 

"0T1“5‘9 

2.440 


,‘T °UT  FRO*-'  2656 


REFERENCE  E.NVl  :-C*  “ENT 

ATTENUATION  = 

-0.7269 

1 

/1“+RP7**2  = 

2.26611 

--  "  -PPESE‘ T  F*:v  I  D0N*‘E.NT 

ATT ET.QAT  ION  - 
/l+RP/**2  = 

-C-.704'3 

2.27707 

REFERENCE  ANTE'  *>A  _ _ ^ 

^APE^T'J^E  admittance  MATCHED  TO  =  C.001C710  +  I  (  0. 0005976  ) 

— - - - — - r — f '  tt- ■  i — 7 - r*  — .^r  -  -  -  ^  . 

w  •  ^wL»»»c  wvj  *r  1  l  w#wwwm«/ww  / 

PRESENT  -  C.C010710  +  I  (  0.0005976  ) 

.  . "'0.00C0479  "h-T-T  -0.00016  =  1  ) 


~®/  */*«?  =-"0.000 1T9 - §  1 +R7*  *7 - =  T;  JO  7F7F" 

®orcLrCr::N  LOSS  =  0.0CC5  ^ISOLATION  =  21.5147 


